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The plant terpenoids limonene (CioHig) and a-bisabolene (C15H24) are hydrocarbon precur- 
sors to a range of industrially relevant chemicals. High-titer microbial synthesis of limonene 
and a-bisabolene could pave the way for advances in in vivo engineering of tailor-made 
hydrocarbons, and production at commercial scale. We have engineered the fast-growing 
unicellular euryhaline cyanobacterium Synechococcus sp. PCC 7002 to produce yields of 
4mgL^^ limonene and 0.6 mgL^^ a-bisabolene through heterologous expression of the 
Mentlia spicata L-limonene synthase or the Abies grandis (E)-a-bisabolene synthase genes, 
respectively. Titers were significantly higher when a dodecane overlay was applied during 
culturing, suggesting either that dodecane traps large quantities of volatile limonene or 
a-bisabolene that would otherwise be lost to evaporation, and/or that continuous prod- 
uct removal in dodecane alleviates product feedback inhibition to promote higher rates of 
synthesis. We also investigate limonene and bisabolene production in the AglgC genetic 
background, where carbon partitioning is redirected at the expense of glycogen biosyn- 
thesis. The Synechococcus sp. PCC 7002 AglgC mutant excreted a suite of overflow 
metabolites (a-ketoisocaproate, pyruvate, a-ketoglutarate, succinate, and acetate) duhng 
nitrogen-deprivation, and also at the onset of stationary growth in nutrient-replete media. 
None of the excreted metabolites, however, appeared to be effectively utilized for terpenoid 
metabolism. Interestingly, we observed a 1.6- to 2.5-fold increase in the extracellular con- 
centration of most excreted organic acids when the AglgC mutant was conferred with 
the ability to produce limonene. Overall, Synechococcus sp. PCC 7002 provides a highly 
promising platform for terpenoid biosynthetic and metabolic engineering efforts. 
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INTRODUCTION 

In recent years, cyanobacteria have been successfully used as a 
platform to generate a range of commercially significant prod- 
ucts, including isoprene (Lindberg et al., 2010; Bentley and Metis, 
2012; Bentley et al, 2014), isobutanol (Atsumi et al, 2009), 2,3- 
butanediol (Oliver et al., 2013, 2014), and ethylene (Lingerer et al., 
20 1 2 ). In all cases, photoassimUated carbon was drawn from native 
metabolic pathways through the heterologous expression of one 
or more enzymes to create new carbon sinks. While demonstrating 
the proof-of-principle, these strategies are unlikely to attain high 
yields of novel end-products as carbon partitioning and regula- 
tion are hardwired into the organism's central metabolic networks. 
For cyanobacteria, one of the main metabolic engineering bottle- 
necks is the ability to redirect carbon from the major glycogen 
and protein sinks to pathways of interest. A careful rewiring of 
central metabolism is required to increase product yield at the 
expense of biomass accumulation, without negatively influencing 
photosynthesis and carbon fixation rates. 

In the context of photoautotrophic biotechnology, the ter- 
penoid biosynthetic pathway is of particular interest, as it produces 



the largest and most diverse array of naturally occurring organic 
compounds (typically of plant origin) (Davies et al., 2014). The 
plant terpenoids limonene (CioHie) and bisabolene (C15H24) are 
recognized as precursors to a range of commercially valuable 
products, with applications in biofuels, bioplastics, pharmaceu- 
tical, nutraceutical, and cosmetic industries (Duetz et al., 2003; 
Peralta-Yahya et al, 2011). In this study, we engineer the fast- 
growing euryhaline cyanobacterium Synechococcus sp. PCC 7002 
for the production of limonene and a-bisabolene through het- 
erologous expression of the Mentha spicata L-limonene synthase 
and Abies grandis (E) -a-bisabolene synthase genes. In cyanobac- 
teria, terpenoids are synthesized via the 2-C-methyl-D-erythritol 
4-phosphate (MEP) pathway, utilizing photosynthetically gener- 
ated pyruvate and glyceraldehyde-3-phosphate (GAP) as the pri- 
mary substrates (Figure 1). Limonene and bisabolene represent 
new carbon sinks that draw from the native terpenoid path- 
way at the levels of geranyl pyrophosphate (GPP) and farnesyl 
pyrophosphate (FPP). 

The two major competing pathways for cyanobacterial ter- 
penoid production are protein and glycogen (the cyanobacterial 
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FIGURE 1 I Terpenoid biosynthesis via tlie 2-C-metliYl-D-erytliritol 
4-pliospliate (IVIEP) pathway. The IVIEP terpenoid biosynthetic pathway 
converts pyruvate and glyceraldehyde-3-phosphate (GAP) feedstock to the 
C5 isomers, isopentenyl pyrophosphate (IPP), and dimethylallyl 
pyrophosphate (DMAPP).The stepwise addition of IPP to DIVIAPP 
generates CIO geranyl pyrophosphate (GPP) and C15 farnesyl 
pyrophosphate (FPP), the respective precursors to monoterpenes and 
sesquiterpenes. In this study, the heterologous expression of the M. 
spicata L-limonene synthase (MsLS) and A. grandis (E)-(!(-bisabolene 
synthase (AgBIS) in Synechococcus sp. PCC 7002 was designed to draw 
on GPP and FPP pools, respectively, generated by the native cyanobacterial 
MEP pathway for the production of the non-native terpenoids L-limonene 
and a-bisabolene. 



Storage carbohydrate) biosynthesis. Recent studies in Synechocystis 
sp. PCC 6803 and Synechococcus sp. PCC 7942 have demonstrated 
that upon simultaneous inactivation of protein and glycogen 
sinks {via nitrogen-deprivation in the AglgC glycogen-deficient 
mutant), a suite of organic acids are excreted due to the result- 
ing metabolic imbalance (Carrieri et al., 2012; Grundel et al, 
2012; Hickman et al., 2013). In our study, we sought to eluci- 
date how such major shifts in flux distribution wUl affect the 
overall metabolic partitioning toward the terpenoid pathway in 
Synechococcus sp. PCC 7002, using the newly engineered limonene 
and a-bisabolene terpenoid sinks as the reporting products. 

MATERIALS AND METHODS 
CLONING AND PLASMID CONSTRUCTION 

The E. coli DHSa strain was used for routine subcloning accord- 
ing to standard protocols. The (— )-4S-limonene synthase gene 
from M. spicata (GenBank accession Q40322) (Colby et al., 1993) 
was optimized for codon-usage in Synechococcus sp. PCC 7002 
(DNA 2.0, USA) without the predicated chloroplast transit pep- 
tide (cTP), and is referred to in this work as MsLS. The first 
168 bases of Q40322 encoding the cTP were replaced with an 
ATG codon for translation initiation in Synechococcus sp. PCC 
7002, retaining residues R58/R59 that are required for catalytic 



activity (Williams et al, 1998). The (£) -a-bisabolene synthase 
gene from A. grandis (GenBank accession 081086) (Bohlmann 
et al, 1998; Trapp and Croteau, 2001) was also optimized for 
codon-usage in Synechococcus sp. PCC 7002 (DNA 2.0, USA), 
and is referred to as AgBIS. Oligonucleotide primers used to 
amplify these codon-optimized genes are as follows (introduced 
restriction sites for cloning are indicated in bold, and lower case 
letters represent start or stop codons): MsLS, MsLS_NcoI_F, 5'- 
CTTCGCTGAACCatgGAACGTCGCTC-3', and MsLS_BamHI_R, 
5'-TGACGGGATCCttaAGCGAAGG-3'; AgBIS, AgBIS_NcoI_F, 
5'-AGATTAATTCCatgGCCGGTGTGA-3';andAgBIS_BamHI_R, 
5'-GCCGGATCCttaTAACGGCAAC-3'. 

A neutral site (NSI) for the integration of transgenes 
was identified in the Synechococcus sp. PCC 7002 genome 
between two open-reading frames encoding hypothetical pro- 
teins (SYNPCC7002_A0935 and SYNPCC7002_A0936), which 
we have termed NSI. This was used as the site of integra- 
tion for the limonene synthase (MsLS) and a-bisabolene syn- 
thase (AgBIS) genes. Oligonucleotide primers designed to amplify 
750 bp up-stream (us) and down-stream (ds) regions flank- 
ing NSI, which were used for double homologous recombina- 
tion in Synechococcus sp. PCC 7002 are as follows (introduced 
restriction sites are indicated in bold): NSIus, 935us_NsiI_F, 
5'-AGTTCACATGCATAAAGTCAA-3', and 935us_EcoRI_R, 5'- 
CGTATAGGAATTCTTACTCAG-3'; NSIds, 936ds_SaII_F, 5'- 
GCATACTGTCGACCTATTTTAT-3'; and 936ds_SphI_R, 5'- 
AGTTGACGCATGCAGAGGTGG-3'. These regions were cloned 
into a pAQl plasmid containing the Synechocystis sp. PCC 6803 
cpcBA (cpc) promoter (Xu et al, 2011). The MsLS and AgBIS 
transgenes, along with a spectinomycin-resistance cassette (SmR), 
were cloned between the NSI-flanking regions immediately down- 
stream of the cpc promoter, generating plasmids pNSI-cpc-MsLS- 
SmR and pNSI-cpc-AgBIS-SmR. To demonstrate that any pheno- 
types arising from the expression of MsLS or AgBIS were not a 
direct result of interference with the NSI region of Synechococcus 
sp. PCC 7002 genomic DNA, we also constructed the control plas- 
mid pNSI-cpc-YFP-SmR, containing a YFP gene in place of the 
MsLIM or AgBIS terpene synthase (TPS) genes. 

CYANOBACTERIAL STRAINS, GROWTH CONDITIONS, AND 
TRANSFORMATION 

Wild type and transformant Synechococcus sp. PCC 7002 strains 
were maintained on solid A+ media (Stevens et al., 1973) supple- 
mented with 8.25 mM Tris-HCl (pH 8.2) and 0.3% sodium thio- 
sulfate. Where appropriate, spectinomycin was used at a concen- 
tration of 50 (ig mL^^ and kanamycin at 100 (ig mL^^ . Liquid cul- 
tures were grown in A+ media supplemented with 8.25 mM Tris- 
HCl (pH 8.2). For nitrogen-deprivation, NaNOa was omitted from 
the growth media and replaced mohmol with NaCl. Nitrogen- 
deplete media is referred to as A+(— N). Cultures were grown in 
250 mL Erlenmeyer flasks on an orbital shaker at 37°C, in an incu- 
bator with an atmosphere of 1% CO2 in air, and constant illumi- 
nation at 250 (xmol photons m^^ sT^ of photosynthetically active 
radiation (PAR). Liquid starter cultures were grown in the presence 
of antibiotics; however, during experimental procedures antibi- 
otics were eliminated. Transformation of Synechococcus sp. PCC 
7002 with the plasmid constructs described above was performed 
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according to established procedures for cyanobacteria (Eaton-Rye, 
2011) to generate the strains LS (containing MsLS) and BIS (con- 
taining AgBIS). Complete chromosomal segregation for the intro- 
duced transgenes was achieved through propagation of multiple 
generations on selective agar and verified by colony-PCR using the 
primers: NSI_us_F, 5'-CTAGCACAAATGAAGCCCGAC-3', and 
NSI_ds_R, 5'-GCAGATATAAGCAACGGTACAG-3' (Figure 2). 
The A.glgC strain was obtained from D. A. Bryant, which was 
generated via the insertional disruption of the glgC open read- 
ing frame (SYNPCC7002_A0095) with a kanamycin-resistance 
cassette (KmR) (Guerra et al., 2013). Oligonucleotide primers 
used to verify glgC disruption in transformant strains are: 
glgC_F, 5'-TCACGTAGTCGGGTTTGATGTC-3', and gIgC_R, 5'- 
CACTAAAGTCCACGACACGACC-3'. The AglgC strain was also 
transformed with the MsLS and AgBIS plasmid constructs to 
generate strains AglgC:LS and Ag/^C:BIS, respectively. 

TERPENOID EXTRACTION AND QUANTIFICATION 

Synechococcus cultures used for terpenoid production assays were 
grown photoautotrophically in 100 mL batches, with an inocula- 
tion density of OD730 nm = 2.0. Where indicated, a 5-mL dodecane 
(analytical grade, Sigma- Aldrich, USA) overlay was applied to the 
culture at f = 0 to trap terpenoids during the 96-h experimen- 
tal time course. At 96 h, cells were harvested by centrifugation to 
separate the biomass from the media. Terpenoid extraction was 
performed on the cell pellet using the chloroform:methanol:water 
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FIGURE 2 I MsLS and AgBIS transgene integration to the 
Synechococcus chromosome. (A)Tlie YFR IVIsLS, and AgBIS transgene 
constructs were introduced to a neutral site (NSI) of genomic DNA 
(between open-reading frames A0935 and A0936) via double homologous 
recombination. Constructs contained the Synechocystis sp. PCC 6803 
cpcBA (cpc) promoter to drive transgene expression, and a 
spectinomycin-resistance selectable marker (SmR). Arrows indicate the 
position of genomic DNA primers used to check for integration and 
complete chromosomal segregation in transformant lines. (B) Verification of 
compete chromosomal segregation for the introduced transgenes in YFR 
LS, and BIS transformants by colony-PCR using the above-mentioned 
primers. Three independent transformant lines of LS (LS.l, LS.2, and LS.3) 
and BIS (BIS.l, BIS.2, and BIS.3) were verified. 



biphasic separation protocol described by Bhgh and Dyer (1959). 
To determine if any accumulated terpenoids had partitioned to the 
media or adhered to the flask walls, the supernatant was returned 
to the original Erlenmeyer flask and a 5-mL dodecane overlay 
applied for 1 h for terpenoid extraction. 

Aliquots of 200 (xL were harvested from the dodecane overlay 
at each experimental time point. These samples, and the chlo- 
roform extracts from the cell pellet, were analyzed directly using 
1 [iL of sample injected into a GC-FID (Agilent 7890A) equipped 
with a DB-5-ms column (Agilent Technologies, Santa Clara, CA, 
USA; 30 m X 0.25 mm; 0.25 [im film thickness) and operated 
using the temperature program: 60-310°C (at 12°Cmin^'), and 
H2 carrier gas at 1.5mLmin^'. Quantification of limonene or 
a-bisabolene from the dodecane or chloroform extracts was per- 
formed using standard curves constructed by GC peak integration 
of serial dilutions of limonene (Sigma- Aldrich, USA) and bis- 
abolene (Alfa Aesar, USA) standards in the appropriate solvent. 
Due to the reverse solvent effect (Grob, 1983; Grob and Schilling, 
1983), limonene was observed as a broad peak on the chro- 
matogram when mixed with dodecane. Dilution of the sample 
with a solvent of lower molecular weight than the solute typi- 
cally eliminates this problem. However, dilution pushed limonene 
below the limits of detection in some samples, and so limonene 
peaks were integrated using the area measured under the broad 
limonene peak. The bisabolene standard contained a number 
of discrete peaks, and -40% of the total combined peak areas 
were attributed to other sesquiterpenes with different retention 
times relative to a-bisabolene. The weight contribution by a- 
bisabolene was adjusted accordingly for the standard curves. Mass 
spectral analysis was conducted using a Varian 3800 GC and Var- 
ian 1200 quadrupole MS/MS equipped with a Rxi-5ms column 
(30 m X 0.25 mm; 0.25 |xm film thickness); temperature program: 
60-310°C (at 12°C min^^); He carrier gas at 1.2 mLmin^^; mass 
spectra, 70 eV, EI mode; ion source temperature 215°C; scan mass 
range, 18-500. 

EXCRETED ORGANIC ACID MEASUREMENTS 

Synechococcus cultures used for excreted metabolite assays were 
grown photoautotrophically in 50 mL batches with an inocula- 
tion density of OD730nm = 2.0. Cells were washed twice with 
A+(— N) media to remove media nitrogen before inoculation in 
either nitrogen-replete or nitrogen-deplete media for experimen- 
tation. Aliquots of cell suspension (ImL) were removed from 
the culture at appropriate time points and centrifuged to pel- 
let cell biomass. The resulting supernatant was filtered through a 
0.45-(xm silicon membrane and the filtrate used directly for high- 
performance liquid chromatography (HPLC) and nuclear mag- 
netic resonance (NMR) spectroscopy. A Surveryor Plus (Thermo 
Scientific, Waltham MA) HPLC equipped with an Aminex fer- 
mentation monitoring column [(150 mm x 7.8 mm); Bio-Rad, 
Hercules, CA, USA] was used for HPLC analyses. 25 \lL of filtrate 
was injected for isocratic elution using 8 mM H2SO4 as the mobile 
phase with a flow rate of 0.5mLmin^^ and a column operating 
temperature of 45°C. A refractive index (RI) detector operating 
at 50°C was used for metabolite identification and quantification. 
Metabolite peaks were integrated and quantified using standard 
curves constructed from serial dilutions of authentic standards 
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(Sigma- Aldrich, USA). Proton NMR spectra were generated with 
a JEOL EGA 500 MHz spectrometer; 128 scans were used with 
water-suppression. 

BIOMASS, CHLOROPHYLL a. CARBOHYDRATE, AND WHOLE CELL 
ABSORPTION SPECTRA MEASUREMENTS 

Cyanobacterial biomass accumulation was measured as dry cell 
weight (DCW), where 2 mL samples of culture were pelleted and 
washed with 8.25 mM Tris-HCl (pH 8.2) to remove any resid- 
ual salt from the growth medium. Pellets were then resuspended 
in 2 mL 8.25 mM Tris (pH 8.2), and dried overnight at 80°C in 
pre-weighed aluminum dishes prior to measuring the DCW. The 
dry weight of 2 mL 8.25 mM Tris (pH 8.2) solution (without 
cells) was subtracted from the cell dry weight to eliminate Tris 
buffer contributions. Chlorophyll a concentrations in cell suspen- 
sions were determined spectrophotometricaUy in 90% methanol 
extracts according to Meeks and Castenholz (1971). Total reduc- 
ing sugars (in glucose equivalents) were quantified by the an throne 
assay, as described by Meuser et al. (201 1 ), as a measure of the car- 
bohydrate content of cells. For whole cell absorption spectra, cell 
suspensions were normalized to an OD730nm of 0.3 and measured 
spectrophotometricaUy (400-800 nm). 

RESULTS 

CONSTRUCTION OF SYNECHOCOCCUS sp. PCC 7002 STRAINS 
CARRYING MsLS AND AgBIS TRANSGENES 

A "NSf region of the Synechococcus sp. PCC 7002 genome was 
identified between two open reading frames encoding hypotheti- 
cal proteins and used as a site for transgene integration (described 
in Materials and Methods). The cpcBA promoter from Synechocys- 
tis sp. PCC 6803 (Xu et al., 2011) was used to drive expression of 
the TPS transgenes. Numerous lines were selected for each new 
strain generated to ensure that any phenotypes arising were a con- 
sistent feature of the specific transformation. Figure 2 shows the 
successful integration of the MsLS and AgBIS constructs at the 
NSI locus, as evidenced by the gel shift of the amplified colony- 
PCR product in transformant strains relative to the wild type. 
Integration of the YFP construct was also confirmed in the YFP 
transformant, which was designed as a control strain to identify 
phenotypes in the LS and BIS strains that may be attributed to 
chromosomal disruption at the NSI locus, rather than expression 
of the TPS transgenes. Importantly, the lack of a wild type-sized 
band in all transformant lines indicated that homoplasmy for the 
introduced transgenes was achieved. 

IDENTIFICATION OF LIMONENE AND (y-BISABOLENE AS 
END-PRODUCTS OF TRANSFORMANT STRAINS 

The use of an organic solvent overlay has proven to be a successful 
method for harvesting terpenoids from microbial cultures (New- 
man et al, 2006; Anthony et al., 2009; Alonso-Gutierrez et al, 2013; 
Bentley et al., 2013). GG-FID analyses of solvent overlays from LS 
(Figure 3A) and BIS (Figure 3B) transformants showed promi- 
nent peaks with similar retention times to commercial standards 
of L-limonene (4.40 min) and bisabolene (9.89 min), respectively, 
which were absent from the wild type. Figure 4 shows the mass 
spectral analyses of putative L-limonene and a-bisabolene peaks 
from transformant solvent extracts in comparison to the corre- 
sponding reference spectra from the NIST Mass Spectral Library 



(upper panels). The solvent also extracted additional products 
common to all strains, including wild type (Figure 3, -10.5 min 
retention time), which may be membrane lipids (accumulated 
levels were below detection limits for mass spectral analysis). A 
comparison of L-limonene and bisabolene commercially available 
standards with NIST reference spectra is also shown in Figure 4 
(lower panels). The major peak (9.89 min) in the bisabolene stan- 
dard (Figure 3B) was identified as a-bisabolene, the peak at 
9.94 min was bisabolol, while the small peaks with retention times 
<9.98 min were aU sesquiterpenes based on a mass to charge (m/z) 
of 204. 

TRANSFORMANT FITNESS AND EFFECT OF DODECANE ON THE RATES 
OF PHOTOAUTOTROPHIC GROWTH 

Photoautotrophic growth rates, as measured by chlorophyll and 
biomass accumulation, were comparable between wild type and 
the YFP control strain (Figure 5A, left panels). This verified 
that NSI is a suitable "neutral" site for transgene integration in 
Synechococcus sp. PCC 7002 under the defined experimental con- 
ditions, because cell growth rate was not adversely affected by 
the genetic manipulation at this locus. Rates of chlorophyll and 
biomass accumulation in the LS and BIS transformants were 
also comparable to that of the wild type (Figure 5, left pan- 
els), confirming that expression of the MsLS and AgBIS trans- 
genes, and subsequent intracellular accumulation of limonene 
and a-bisabolene, did not adversely affect cell fitness under our 
experimental conditions. 

Dodecane was chosen as the solvent of choice for longer- 
term terpenoid extraction experiments due to its relatively low 
volatility, which enabled continuous extraction over multiple days 
(evaporative loss was determined as negligible over the five day 
cultivation period, data not shown). The Cu chain length also 
allowed chromatographic separation from Gio limonene and G15 
bisabolene. The rate of limonene and a-bisabolene biosynthesis 
during the growth phase was measured over a 96-h time course 
using a 5% (v/v) dodecane overlay. Notably, the presence of the 
dodecane overlay did not greatly affect chlorophyll content or bio- 
mass yield in any of the strains (Figure 5, right panels), indicating 
that atmospheric CO2 was able to efficiently diffuse though the 
dodecane layer and dissolve in the liquid media at concentrations 
that were not limiting for photoautotrophic growth under these 
conditions. 

LIMONENE AND BISABOLENE BIOSYNTHETIC RATES AND CELLULAR 
LOCALIZATION 

Three independent lines of the LS transformant consistently pro- 
duced yields over 4 mg limonene L culture^ ^ , with the highest rate 
of 50|J-gL culture^^h^^ recorded over the exponential growth 
phase when cells were actively dividing (Figure 6A). As cultures 
reached higher cell densities, light limitation due to cell shading 
began to limit photosynthetic growth (Figure 5A), which likely 
translated into reduced rates of limonene biosynthesis. Assuming 
that carbon comprised -50% of the total biomass, the carbon 
partitioning to limonene was estimated at 0.3%. Yields of a- 
bisabolene were considerably lower, where average yields in BIS 
transformant lines ranged between 0.5 and 0.7 mg a-bisabolene L 
culture^' (Figure 6B). The highest a-bisabolene production rate 
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FIGURE 3 I GC-FID analyses of solvent overlay-extracted limonene and 
bisabolene from transformant strains. (A) Hexane-overlay extraction from 
tlie LS transformant (LS.l line as a representative) showed accumulation of a 
product that was absent in the wild type, and with a similar retention time 
(4.40 min) to a commercial L-limonene standard. (B) Dodecane overlay 
extraction from the BIS transformant (BIS.l line as a representative) revealed 



a peak with a similar retention time (9.89 min) to a bisabolene standard, which 
is absent in the wild type. For these experiments, dodecane was used for 
bisabolene extraction because of the extended period of time required to 
accumulate detectable levels of bisabolene, whereas hexane was applied for 
better limonene peak resolution due to the reverse solvent effect observed 
with dodecane. 



was 6 |xg L culture^ over the active growth phase, equivalent 
to -0.06% of assimilated carbon partitioning as a-bisabolene. 

A methanohchloroform extraction was performed on cells har- 
vested at the termination of the 96-h time course to assess whether 
any limonene or a-bisabolene remained in the cells after the 
dodecane extraction procedure. On average, 1 1 (ig limonene L 
culture^' was extracted from LS cells (representing 0.25% of the 
total limonene yield), and 40 |ig a-bisabolene L culture^' from 
BIS cells (8% of the total bisabolene yield) (Figure 6C). Therefore, 
almost 2.5-fold (mohmol) more a-bisabolene remained inside 
cells after dodecane extraction, compared to limonene. Relative 
to the overall terpenoid yields from LS and BIS transformants, a 
much greater fraction of a-bisabolene was retained intraceUularly 
compared to limonene. This may be explained, in part, by the 
mass difference between these two compounds, with the smaller 
limonene molecule having a greater rate of diffusion across the 
cell membrane. 

Surprisingly, similar amounts of limonene and a-bisabolene 
were extracted from harvested cells regardless of the presence or 
absence of the dodecane overlay (Figure 6C). Additional dodecane 
extractions of the spent media performed after cell harvesting did 
not reveal any terpenoids that partitioned to the media. The higher 



overall product yields displayed by cultures overlayed with dode- 
cane maybe due to: (i) the trapping and solubilization of limonene 
and a-bisabolene, which prevented evaporative loss, and/ or (ii) the 
sequestration of limonene and a-bisabolene from the cells, which 
alleviated a negative feedback inhibition mechanism. 

PHYSIOLOGICAL CHARACTERIZATION OF SYNECHOCOCCUS sp. PCC 
7002 AglgC 

At a light intensity of 250 [xmol photons rcT^ s^^ PAR, the A.glgC 
strain displayed slightly impaired photoautotrophic growth in 
nutrient-replete media, as measured by ODysonm and chlorophyll 
content (Figure 7 A, left panels). In contrast to the wild type, A.glgC 
did not grow photoautotrophicaUy in nitrogen-deplete medium 
(Figure 7 A, right panels) due to its inability to degrade the light- 
harvesting phycobilisomes as a source of nitrogen (Carrieri et al, 
2012; Grundel et al., 2012; Guerra et al, 2013; Hickman et al, 
2013). Accordingly, the absorbance spectrum of AglgC cells in 
nitrogen-deplete media shows the presence of phycobilin, which 
is absent in the wild type after nitrogen-deprivation (Figure 7B). 
As a result, the AglgC culture retained the blue-green hue under 
nitrogen-deprivation, while the wild type appeared more yellow- 
green due to the unmasking of chlorophyll a as the blue-pigmented 
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FIGURE 4 I Mass spectra analyses of limonene and a-bisabolene. 

(A) Comparison of the L-limonene reference spectra from the NIST library 
with the spectra of the putative L-limonene peak in the LS.l line (upper panel) 



and the authentic L-limonene standard (lower panel). (B) a-Bisabolene NIST 
reference spectra alignment with the putative a-bisabolene peak in the BIS.l 
line (upper panel) and a commercial bisabolene standard (lower panel). 



phycobiliproteins were degraded (Figure 7C). We confirmed that 
the reducing carbohydrate content of AglgC was diminished 
compared to the wild type, particularly during stationary phase 
at 48 h in nutrient-replete media, and also during the first 24 h 
after the onset of nitrogen-deprivation (Figure 7D). These data 
confirmed that AglgC is defective in glycogen biosynthesis, and 
that our experimental conditions replicated the phenotypes previ- 
ously observed in AglgC mutants of other cyanobacterial species 
(Carrieri et al, 2012; Grundel et al., 2012; Guerra et al., 2013; 
Hickman et al., 2013), and in Synechococcus sp. PCC 7002 (Guerra 
etal, 2013). 

ORGANIC ACID EXCRETION IN AglgC AND ITS EFFECT ON TERPENOID 
YIELDS 

Consistent with studies in Synechocystis sp. PCC 6803 (Carri- 
eri et al., 2012; Grundel et al., 2012) and Synechococcus elonga- 
tus PCC 7942 (Hickman et al., 2013), we identified pyruvate, 
a-ketoglutarate, and succinate in the spent media of nitrogen- 
deprived Synechococcus sp. PCC 7002 AglgC cultures using HPLC 
analysis (Figure 8A, left panels), through a comparison of elu- 
tion times with known standards (Figure 8B). However, we also 
detected a-ketoisocaproate and acetate within the same concen- 
tration range, metabolites which were not previously reported in 
the AglgC strains of the other cyanobacterial species. Acetate is a 
common fermentation product putatively derived from acetyl- 
CoA in cyanobacterial metabolism (Xu et al., 2013), while a- 
ketoisocaproate is the immediate precursor to leucine. This high- 
lights subtle metabolic distinctions in the metabolism of the AglgC 
mutants between the model cyanobacterial species. An additional 
unique observation was the excretion of the same organic acids by 
AglgC in nutrient-replete media at the onset of stationary phase, 
when wild type cells would normally begin to increase glycogen 



stores for energy reserves (Figure 8A, right panels). However, on a 
per-cell basis, these concentrations were much lower relative to the 
concentrations of those excreted during nitrogen starvation. NMR 
was used as a second, independent method to verify the identity 
of these metabolites (Figure 8C). Pyruvate, succinate, and a- 
ketoisocaproate accumulated at concentrations of -250 (j,M, while 
a-ketoglutarate accumulated at -200 ^lM, and acetate -500 [xM 
under nitrogen starvation. 

To assess whether the metabolic flux imbalance in AglgC 
redirected metabolites toward terpenoid biosynthesis, we grew 
the A^/gC:LS and Ag/^C:BIS transformants under conditions of 
nitrogen starvation. The total yields of limonene and a-bisabolene 
were extremely low relative to the culture volume, which was not 
surprising given the absence of growth during nitrogen starvation 
(data not shown). However, terpenoid yield relative to biomass 
was relatively constant, regardless of whether cultures were grown 
in the presence or absence of nitrogen. In fact, throughout this 
study, terpenoid yield as a function of biomass remained relatively 
constant among all combinations of genetic background {±glgC), 
nutrient availability (initrogen), and growth phase (exponential 
vs. stationary). This indicates that terpenoid biosynthesis is cou- 
pled to cell growth and that flux through the terpenoid pathway is 
tightly regulated, even in the face of metabolic perturbations that 
increase feedstock availability. 

Interestingly, however, we observed a higher extracellular con- 
centration of excreted organic acids in the limonene-producing 
AglgC:LS strain, relative to the AglgC and AglgC'.BlS strains 
(Figure 9). Under nitrogen starvation, Ag/^C:LS accumulated 
-2.5-fold more a-ketoglutarate, and almost twofold more a- 
ketoisocaproate, pyruvate, and succinate in the spent media com- 
pared with the AglgC strain, despite no detectable photoau- 
totrophic growth in either strains. It is unclear whether the higher 
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FIGURE 5 I PPiotoautotrophic growth of transformant lines in the 
presence and absence of a dodecane overlay. (A) Chloropliyll and 
biomass accumulation las measured by dry cell weight (DCW)] over 96 h in 
the YFP control strain, and three independent LS lines (LS.1 , LS.2, and LS.3) 
in comparison to wild type. Left panels show growth in the absence of a 
solvent overlay, while right panels display growth in the presence of a 
dodecane overlay. (B) Comparison of chlorophyll and biomass accumulation 
(DCW) between wild type and three independent BIS transformant lines 
(BIS.l , BIS. 2, and BIS. 3), in the absence (left panels) and presence (right 
panels) of a dodecane overlay. Error bars represent standard deviation from 
at least three biological replicates. 



concentrations in the AglgC:LS strain are due to increased metab- 
olism toward these end-products, or if a greater proportion of the 
organic acids are excreted from the cell, but provides an interesting 
line of future investigation. 

DISCUSSION 

We have demonstrated for the first time, the synthesis of limonene 
in Synechococcus sp. PCC 7002, and the first synthesis of a- 
bisabolene in a photosynthetic microorganism. The choice of 
M. spicata (— )-4S-limonene synthase and A. grandis {E)-a- 
bisabolene synthase as representative monoterpene and sesquiter- 
pene synthases was based on proven successes in heterotrophic 
microbial production hosts (Carter et al., 2003; Peralta-Yahya et al., 
2011; Alonso-Gutierrez et al., 2013; Ozaydin et al, 2013). 

Importantly, (— )-4S-limonene synthase and (£)-()( -bisabolene 
synthase are known to have high product specificity (Alonso 
et al, 1992; Rajaonarivony et al, 1992; Bohlmann et al, 1998), 
which is unusual for TPSs. Most TPSs generate multiple prod- 
ucts from the prenyl diphosphate substrate as a result of the 
highly unstable carbocationic intermediates generated during the 
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FIGURE 6 I Rates of limonene and cc-bisabolene production in 
transformant lines. (A) Limonene accumulation in three Independent LS 
transformant lines (LS.l, LS.2, and LS.3). Yields were quantified by GC-FID 
analysis of limonene extracted by a dodecane overlay during the 96-h 
growth phase (Figure 5A, -l-dodecane). (B) a-Bisabolene extracted from 
three independent BIS lines (BIS.l, BIS.2, and BIS.3) over the 96-h growth 
phase (Figure 5B, +dodecane) using a dodecane overlay. (C) Yields of 
limonene and a-bisabolene extracted from cells harvested after 96 h of 
photoautotrophic growth (Figure 5) in the presence (solid bars) or absence 
(hatched bars) of a dodecane overlay. Terpenoids were extracted from cells 
using a methanohchloroform extraction method, and quantified by GC-FID. 
Error bars represent standard deviation from at least three biological 
replicates. 



electrophilic reaction. Accordingly, these enzymes also produced 
pure L-limonene and a-bisabolene products in Synechococcus sp. 
PCC 7002; however, the yield of L-limonene was much greater 
than that of a-bisabolene. This could be explained by the supe- 
rior enzyme kinetics of the M. spicata L-limonene synthase, which 
has a (GPP) of 1.8 (i-M and k^^ of 0.3 s"' (Alonso et al, 
1992; Rajaonarivony et al, 1992), compared to the A. grandis 
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FIGURE 7 I Physiologv of Synechococcus sp. PCC 7002 AglgC. 

(A) . Comparison of photoautotrophic growth rates in wild type and AglgC 
cultures in nitrogen-replete and nitrogen-deplete media, as measured by 
changes of optical density (730 nm) and chlorophyll content over 48 h. 

(B) Whole cell absorption spectra of wild type and AglgC cells after 48 h 
growth in nitrogen-replete and nitrogen-deplete media. (C) Differences in 
pigmentation of wild type and AglgC cultures after 48 h growth in the 
presence or absence of nitrogen in the media. (D) Carbohydrate content 
(reported as glucose equivalents) of wild type and AglgC cells at the 
different stages of growth (0, 24, and 48 h) shown in (A). Error bars 
represent standard deviation from three biological replicates, and are 
hidden beneath the marker if not apparent. 



(E)-ci(-bisabolene synthase, with a high K^n (FPP) of 49.5 (iM 
and very low kcut of 0.11 min^^ (McAndrew et al., 2011). Sub- 
strate availability is also likely to play a significant role because 
monoterpenes and sesquiterpenes are absent in most cyanobac- 
teria, and so GPP and FPP substrates may not accumulate to the 
levels seen in plants. However, GPP must accumulate in significant 
quantities as the precursor to GGPP, from which longer chained 
cyanobacterial terpenoids are synthesized (e.g., carotenoids), and 
so may also explain the greater L-limonene yield relative to 
a-bisabolene. 

The cyanobacterium Synechocystis sp. PCC 6803 has previ- 
ously been engineered to produce fi-phellandrene monoterpenes 
(Bentley et al, 2013) and fi-caryophyllene sesquiterpenes 



(Reinsvold et al., 2011). The yield of f5-phellandrene was 
133 (iggDCW^^, as measured by a non-continuous heptane- 
overlay extraction method. Intracellular titers of fi-caryophyllene 
were ~4|xggDCW^', as measured by a methanolxholoroform 
extraction method. However, these yields may have been under- 
estimated due to the volatile nature of the terpenoid products. 
Evaporation was attributed to major losses of the amorphadi- 
ene sesquiterpene from E. coli liquid cultures, and a dodecane 
overlay was demonstrated to reduce the loss (Newman et al, 
2006). Recently, the nitrogen-fixing cyanobacterium Anabaena 
sp. PCC 7120 was engineered for limonene production via het- 
erologous expression of the Picea sitchensis limonene synthase 
(Halfmann et al., 2014). The total yield of limonene (114 jigL^') 
was found to partition as a volatile gas in the culture headspace, 
and titers were improved -ninefold upon overexpression of three 
rate-limiting enzymes involved in terpenoid metabolism. Here, 
we report the yields of 1.7mggDCW^^ (4mgL^^) of limonene 
and0.3mggDCW^' (0.6 mgL^') of a-bisabolene in Sy«ec/iococ- 
cus sp. PCC 7002, upon expression of the single TPS genes, and 
using a continuous dodecane overlay extraction method that may 
reduce product loss due to evaporation. 

It is also conceivable that dodecane may actively sequester ter- 
penoids from the cell, thereby relieving any negative feedback 
pressures exerted by the products thus promoting higher end- 
product yields. Indeed, a 20% (v/v) dodecane overlay was found 
to alleviate growth inhibition of Synechocystis sp. PCC 6803 that 
was induced by the addition of commercial farnesene at con- 
centrations as low as 0.04% (v/v) (Hellier et al., 2013). Further 
support for this theory was demonstrated in E. coli, where het- 
erologous expression of an export pump in a limonene-producing 
strain increased limonene yield by -60% (Dunlop et al, 2011). 
The microbial toxicity of bisabolene appears to be very low, with 
concentrations of commercial bisabolene up to 20% (v/v) exerting 
no effect on cell growth in E. coli and yeast (Peralta-Yahya et al., 
201 1). Although toxicity information is useful in this respect, the 
product feedback inhibition may not always manifest as growth 
inhibition. To fully understand the increased limonene and a- 
bisabolene yields in the presence of dodecane, investigations are 
required: (i) to probe terpenoid (or other) pathway inhibition by 
these end-products as well as (ii) to quantify the biologically pro- 
duced limonene and bisabolene lost via evaporation in the absence 
of dodecane. It is important to establish whether the product nat- 
urally separates from biomass and how the physical properties of 
the product (mass, boiling point) influence this, or whether the 
product exerts feedback inhibition so that engineering efforts may 
be targeted toward this. 

Targeted metabolic flux redistribution among the major carbon 
sinks is the key to producing novel end-products from engineered 
microbes at commercially available quantities. The goal remains 
to increase product yield at the expense of biomass accumulation 
(cell growth and division), without negatively impacting metab- 
olisms associated with the assimilation and direction of carbon 
to the product. Although disruption of glycogen biosynthesis in 
Synechococcus sp. PCC 7002 produced an oversupply of central 
metabolites (a-ketoisocaproate, a-ketoglutarate, pyruvate, succi- 
nate, and acetate) (Figure 10), there was no apparent increase in 
flux though the MEP terpenoid pathway as conveyed by yields 
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FIGURE 8 I High-performance liquid chromatography and NMR 
identification of secreted organic acids from AglgC. 

(A) HPLC-generated diromatograms showing tlie presence of 
secreted metabolites in AglgC, which are absent from the wild type 
culture media, when grown for 48 h under nitrogen-deplete (left 
panels) and nitrogen-replete (right panels) conditions. (B)The 



retention times of the metabolites is shown relative to a mixture of 
4mM a-ketoglutarate, pyruvate, succinate, acetate, and 
a-ketoisocaproate standards. (C) Proton NMR spectrum of the 
secreted metabolites in AglgC, showing chemical shifts and peak 
splitting patterns that correspond to a-ketoglutarate, pyruvate, 
succinate, acetate, and a-ketoisocaproate. 
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FIGURE 9 I Accumulation of excreted organic acids in the spend media of cells grown under nitrogen starvation for 48 h. Concentrations of organic 
acids were measured by HPLC analysis and presented relative to cell optical density at 730 nm. Error bars represent standard deviation from three biological 
replicates. 



of limonene and a-bisabolene reporter products. Furthermore, 
such a severe perturbation of central carbon metabolism had a 
detrimental effect on the organism fitness as manifested by a 
complete growth inhibition. To that end, introduction of alter- 
native routes for carbon partitioning have proven successful in 



cyanobacteria as carbon fixation rates are stimulated by providing 
alternative sinks for photosynthetically derived NADPH and ATP 
(Ducat et al, 2012; Li et al, 2014). While the MEP terpenoid 
pathway utilizes both NADPH and ATP as cofactors, growth 
and pathway optimization are required so that photosynthesis 
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FIGURE 10 I Metabolic redistribution in Synechococcus sp. PCC 
7002 upon simultaneous blocks to protein and carbohydrate sinl<s. 

Protein (amino acids), carboliydrates (mainly glycogen), and lipids 
represent the major metabolic sinks in Synechococcus sp. PCC 7002. 
When metabolic flux to glycogen and amino acids were simultaneously 
blocked, through inactivation of AglgC and elimination of nitrogen from 



the culture medium (red crosses), a suite of organic acids were secreted 
from the cell as overflow metabolites from central metabolism (red 
boxes). There is potential to redirect these spill-over intermediates 
toward biotechnologically relevant pathways to increase yields of 
engineered products, such as limonene and a-bisabolene via the 
terpenoid biosynthetic pathway. 



and terpenoid biosynthesis remain active despite major metabolic 
redistributions. 
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